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Spin-dependent transport through two coupled single-level quantum dots weakly connected to fer- 
romagnetic leads with collinear magnetizations is considered theoretically. Transport characteristics, 
including the current, linear and nonlinear conductance, and tunnel magnetoresistance are calcu- 
lated using the real-time diagrammatic technique in the parallel, serial, and intermediate geometries. 
The effects due to virtual tunneling processes between the two dots via the leads, associated with 
off-diagonal coupling matrix elements, are also considered. Negative differential conductance and 
negative tunnel magnetoresistance have been found in the case of serial and intermediate geometries, 
while no such behavior has been observed for double quantum dots coupled in parallel. It is also 
shown that transport characteristics strongly depend on the magnitude of the off-diagonal coupling 
matrix elements. 
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I. INTRODUCTION 

Transport properties of double quantum dots (DQDs) 
have recently attracted considerable attention from both 
experimental and theoretical sides 10 1 1 1 This 
is mainly due to the fact that DQDs are one of the 
simplest model systems that mimic behavior of real 
molecules, and are thus frequently referred to as ar- 
tificial molecules. Moreover, double quantum dots 
are considered to play an important role in quan- 
tum computation^ 3 , and spintronics i 12 ' 13 ' 14 They ex- 
hibit a variety of different phenomena, including the 
Pauli spin blockade^ formation of molecular states^ 
spin filtering effects^ or various interference effects, 
such as Fano or Dicke resonances! 15 ' 16 In addition, very 
recently it was shown theoretically that, when cou- 
pled to ferromagnetic leads, double quantum dots dis- 
play a considerable tunnel magnetoresistance (TMR) 
and spin accumulation effects i 15 ' 17 ' 18 ' 19 Spin-dependent 
transport properties of quantum dots have been so 
far mainly addressed in the case of single quantum 
dots,2&2ii22i22£ii2£ This field is already rather well estab- 
lished and transport through single quantum dots cou- 
pled to ferromagnetic leads has been extensively stud- 
ied experimentally^ 6 -^^ 8 -^^! 3 ^^^ 3 - On the other hand, 
theoretical investigations of spin effects in multi-dot 
structures are in relatively initial stage, and so is exper- 
imental implementation of DQDs coupled to ferromag- 
netic leads, which still remains a challenge. 

In this paper we consider the spin-dependent trans- 
port properties of double quantum dots focusing on the 
weak coupling regime. Conductance of the system is 
then determined mainly by discreteness of the dots' en- 
ergy levels and Coulomb correlations, which may lead to 
the Coulomb blockade effect and step-like current-voltage 
characteristics. To calculate the transport characteristics 



in the linear and nonlinear response regimes, we employ 
the real-time diagrammatic technique^ This technique 
allows us to take into account the interference effects re- 
sulting from virtual processes between the two quantum 
dots and the leads as well as renormalization of the dot 
levels. In particular, taking into account the first-order 
self-energies, we calculate the current, conductance and 
tunnel magnetoresistance for various geometries of the 
double quantum dot system. In particular, we analyze 
the transport characteristics in the case of DQDs con- 
nected in series, in parallel, as well as for some inter- 
mediate geometries. We show that the interference ef- 
fects associated with off-diagonal matrix elements of the 
self-energy can significantly influence transport proper- 
ties of the system for parallel and intermediate geome- 
tries. When the quantum dots are coupled in series or 
are in an intermediate geometry, we find negative differ- 
ential conductance (NDC) and negative TMR in some 
transport regimes. These features appear in transport 
through DQD systems. However, they were not found 
in transport through a single quantum dot connected to 
ferromagnetic leads in the corresponding range of param- 
eters. Furthermore, we also analyze the dependence of 
transport properties on the magnitude of the off-diagonal 
matrix elements. Finally, we note that in previous theo- 
retical considerations ; 9 ' 10 ' 18 virtual first-order tunneling 
processes have not been taken into account as they be- 
come relevant for parallel and intermediate geometries of 
double quantum dot systems. 

The paper is organized as follows. In Sec. 2 we de- 
scribe the model of a double quantum dot and outline 
the method used in calculations. Numerical results on 
the current, conductance and tunnel magnetoresistance 
for DQDs coupled in serial, in parallel and for interme- 
diate geometries are presented and discussed in Sec. 3. 
The main focus here is on negative differential conduc- 
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tance and negative tunnel magnetoresistance. Summary 
and final conclusions are given in Sec. 4. 

II. THEORETICAL DESCRIPTION 
A. Model 

We consider two coupled single-level quantum dots 
connected to ferromagnetic leads, as shown schematically 
in Fig. [1] The magnetizations of the leads are assumed to 
be collinear, and the system can be either in the parallel 
or antiparallel magnetic configuration. The system can 
be switched from one configuration to the other by apply- 
ing a weak external magnetic field and sweeping through 
the hysteresis loop, provided the leads have different co- 
ercive fields. The Hamiltonian of the double quantum 
dot system is generally given by 

H = i^loads + HoQD + ^tunnel) (1) 

where the first term, i?i ca ds, describes the left (L) 
and right (R) electrodes in the non-interacting quasi- 
particle approximation, i/i oa ds = Hi, + Hr , with Hp = 
J2wa £ l3kacl ka c l3]l<J (for (3 = L,R). Here, Cg kff (c^) is 
the creation (annihilation) operator of an electron with 
the wave vector k and spin a in the lead /3, whereas 
f/3kcr denotes the corresponding single-particle energy. 
The second term of the Hamiltonian describes the double 
quantum dot and is given by, 

ia i 

+ Uo(nn +n u )(n 2T +n 2 J, (2) 

where a = — <j, ni a = d\ a di a is the particle number op- 
erator for spin a in the dot i (i — 1,2), d\ a {di a ) is 
the respective creation (annihilation) operator, and £i a 
denotes the spin-dependent discrete energy level of the 
i-th dot. Double occupation of the dot i is associated 
with the intradot charging energy Ui, whereas simulta- 
neous occupation of both dots with one electron per dot 
costs the interdot charging energy Uq. In the follow- 
ing we assume U\ = U2 = U, and note that Uo < U 
for typical lateral double dot structures.— We further pa- 
rameterize the quantum dot energy levels by their av- 
erage energy, E a = (e\ c + £2<r)/2 and their difference, 
AE = E\ a — £2<r, respectively, so that E\ a = E a + AE/2 
and £2cr = E a — AE/2. Here we assumed that the dots 
have equal g-factors, then AE is independent of spin even 
in the presence of external magnetic field. Furthermore, 
we also assume E^ = E± = E, if not stated otherwise. 
Apart from this, we assume that the bare energy levels 
of the dots are independent of the applied transport volt- 
age. This can be achieved for instance with suitable gate 
voltages. 

The last term of the Hamiltonian, Eq. (JTJ) , consists of 
two different terms, -f/tunnei = + H t . The first one de- 
scribes the spin-dependent tunneling processes between 




FIG. 1: (Color online) Schematic picture of the DQD sys- 
tem coupled to ferromagnetic leads. The parameter 7^ (for 
j3 — L,R,a =t> 4) describes here a spin-dependent dot-lead 
couplings, whereas a takes into account difference in the cou- 
pling of a given electrode to the two dots (a £ (0, 1)). In 
particular, for a = double quantum dots are in the serial 
geometry, while for a = 1 the system is in the parallel geom- 
etry. 

the quantum dots and external magnetic leads and is 
given by 

^v = EE (Vfuc^di. + V^dlcp^) , (3) 

/3i ker 

where Vp, are the relevant tunneling matrix elements be- 
tween the lead (i and dot i. The second term of H tunne i 
corresponds to hopping between the two quantum dots 
and reads 

ff t = -t^(4A+4A). (4) 

a 

The inter-dot hopping parameter t is assumed to be real 
and independent of the electron spin orientation. We 
also assume that all tunneling processes in the system 
are spin-conserving. 

Due to the coupling to external leads, the dot levels ac- 
quire finite widths. The dot-leads coupling is described 
generally by Fg- = 2npgVpiVg '., where is the den- 
sity of states of the lead (3 for spin a, a = +(— ) for the 
majority (minority) spin electrons, r^ - describes the 
spin-dependent hybridization between the local dot lev- 
els (i, j = 1,2) and the leads, and is directly related to the 
coupling strength between the dots and leads. In princi- 
ple, the coupling parameters may be energy-dependent. 
However, for transport regimes considered in this paper, 
it is well justified to assume that the couplings are con- 
stant within the electron bandJ^ For the considered sys- 
tem, the coupling parameters can be conveniently written 
in a matrix form as 

r a - f r ^ n r <3 12 I (^\ 

\ 1 021 1 /322 / 

where the tunneling matrix elements Vpi are assumed 
to be real and constant, while T^ 12 = T^ 21 = 

1p(^pii^P22) 1 ^ 2 ■ The off-diagonal matrix elements of 
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Tp are associated with various interference effects result- 
ing from virtual first-order tunneling processes between 
the two quantum dots through the states in the leads. 
These off-diagonal matrix elements may be significantly 
reduced in comparison to diagonal matrix elements TZ U . 
Furthermore, for complete destructive interference these 
matrix elements may be totally suppressed. To take the 
interference effects into account, we have introduced the 
parameters <7l and or .. (For calculation of the parameters 
ql and qr see Ref. j35|.) We further assumed that qp are 
real positive numbers and fulfill the condition qp < 1. 
Moreover, by introducing the spin polarization of lead /?, 
P/3 = (pp ~~ Pp)/(pp + Pp )■> the coupling constants in the 
parallel configuration can be simply written as 

for the coupling to the left electrode and 

r TU) _ a ( a o Rv ^ \ m 

for coupling to the right lead. In the above expressions 
7£ = (1 ±p L )T L and 7^ = (I ±p R )T B . Here, we as- 
sume that the couplings are symmetric, Tl = Tr = T/2, 
and a takes into account the difference in the coupling 
of a given electrode to the dots, see Fig. [TJ In princi- 
ple, the parameter a could be different for the left and 
right leads. However, we assume here that the system 
is symmetric, as shown in Fig.l. In the above formulas 
we have also assumed that in the parallel configuration 
the spin-| (spin- J,) electrons belong to the majority (mi- 
nority) electron bands of the leads. In the antiparallcl 
configuration the couplings are given by Eqs. ([5]) and (J7]) 
with <-> — pr. By varying the parameter a, one can 
change the geometry of the system from serial one for 
a = to the parallel geometry for a = I. For intermedi- 
ate values of a the system is in an intermediate geometry, 
where each of the two dots is coupled to both leads, see 
Fig. [TJ It is worth noting that investigating the effect 
of geometry on transport properties by varying the pa- 
rameter a is certainly relevant from experimental point 
of view. 



B. Method 

In order to determine the transport properties of 
the system we employ the real-time diagrammatic 
technique^ This technique is based on the perturbation 
expansion of the reduced density matrix and the relevant 
operators with respect to the coupling strength T. Wc 
calculate the reduced density matrix p for the double- 
dot system by integrating out the electronic degrees of 
freedom in the leads. The time evolution of p is then 
described by the Liouville equation of the for m 10 ' 34 

ih^-p = [H DQD + H t , p] + T,p . (8) 



The commutator represents the internal dynamics in the 
double dot, which mainly depends on the level separa- 
tion AE and the interdot coupling t. The second part 
of Eq. ([5]) accounts for the tunnel coupling between the 
double dot and external reservoirs. The complex ten- 
sor £ is associated with tunneling processes and tunnel- 
induced energy renormalization of the dot levels. The 
elements of the reduced density matrix are defined as 
Pj* 1 = (xi|/3|X2), where xi an d X2 denote the eigenstates 
of the DQD system. Then, the Liouville equation for 
stationary reduced density matrix can be written in the 
formal 

= (xi|[Hdqd + H t ,p] X2 ) + J2 Klxl P xi ■ ( 9 ) 

Xix'2 

Here, E XlX , 1 denotes the self-energy corresponding to evo- 

X2X2 

lution forward in time from state \x'\) to state \xi) and 
then backward in time from state \x2) to state |x' 2 ). 
The diagonal elements of the reduced density matrix, P£ 
(xi = X2 = x)i correspond to probability of finding the 
DQD system in the state |x). To solve Eq.9 for density 
matrix elements one usually performs a perturbation ex- 
pansion with respect to the dot-lead coupling strength 
r. Then, each term of the expansion can be visualized 
graphically as a diagram defined on the Keldysh contour, 
where the vertices are connect by lines corresponding to 
tunneling processes. The self-energies in respective order 
of expansion can be calculated using the diagrammatic 
rules M- 

In our considerations we take into account the limit 
of weak tunnel coupling between the dots. For serial 
geometry of the double dot system, tunneling between 
the two dots becomes then a bottle-neck for transport 
and may considerably alter the spin-dependent transport 
characteristics of the system, as presented in the next 
section. This is contrary to previous theoretical studies 
of spin-dependent transport in DQDs,— where the hop- 
ping between the two dots was relatively large and trans- 
port took place through highly hybridized molecular-like 
states of the system. To determine transport character- 
istics in the weak coupling regime we perform systematic 
perturbation expansion with respect to the coupling pa- 
rameter r. Furthermore, we assume T <ti ksT. The 
current is then mediated mainly by first-order (sequen- 
tial) tunneling processes, while the higher-order coherent 
tunneling events play a minor role, and it is justifiable 
to neglect them.— ^ We thus investigate the basic trans- 
port properties using the sequential tunneling approxi- 
mation, i.e. we need to determine only the lowest-order 
self-energies which involve one tunneling line. Some ex- 
amples of first-order diagrams relevant for the present 
calculation are shown in the Appendix. 

After calculating the density matrix elements from 
Eq. ([5]), one can determine the sequential current flow- 
ing through the double dot system from the following 
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formula 

X1X2 

X1X2 

where S 1 , \ } denotes the first-order self-energy in which 

X2X2 bJ 

one vertex was substituted by a vertex representing 
the current operator, I = (I R — /l)/2, with Ip = 

-i(e/h) ~£i Eka(^i4k<T d ^ - V^d\ a c^ a ). 

In our analysis we assume that the intra-dot charging 
energy U is relatively large for both dots, much larger 
than the interdot Coulomb correlation energy Uq. Thus, 
only the zero, one and two-particle DQD states are rele- 
vant for transport. Furthermore, in the case of U 3> Uq, 
the occupation probability of the states with two elec- 
trons in the same dot is also vanishingly small. However, 
these states are taken into account as intermediate (vir- 
tual) states in the calculation, providing a natural high- 
energy cutoff. Finally, we would like to emphasize that 
by using the real-time diagrammatic technique we are 
able to take into account the first-order virtual tunnel- 
ing processes between the two dots in a fully systematic 
way.i^i 



III. NUMERICAL RESULTS 

In this section we present and discuss numerical results 
on the charge current, differential conductance and tun- 
nel magnetoresistance of double quantum dots coupled in 
serial, parallel, and intermediate geometries. The TMR 
effect results generally from spin-dependent dot-lead tun- 
neling processes, which in turn leads to the dependence of 
transport characteristics on magnetic configuration of the 
system. The TMR is quantitatively described by the ra- 
tio TMR = (Ip — Iap)/Iap, where Ip and Iap denote the 
currents flowing through the system in the parallel and 
antiparallcl magnetic configurations, respectively 36 i 37 

In the numerical analysis we assume spin degener- 
ate dot levels, Ei a = Ei (for i = 1,2 and a =T,4)- 
We also assume that external electrodes arc made of 
the same ferromagnetic material, p^ = pr = p, and 
that the system is symmetrically coupled to the leads, 
Tl = Tr = r/2. The parameters (?l and gR can be gen- 
erally different. However, we assume that they are real 
and equal, q^ — qpi = q. We also set the intra- and 
inter-dot Coulomb parameters to be: U = lOOfcsT and 
Uq = 20fcsT, respectively. Finally, we assume spin po- 
larization p = 0.4, which is typical of 3d ferromagnetic 
metals^ The inter-dot hopping parameter t is assumed 
to be: t = 0.25r with T = 5/xeV. These are typical exper- 
imental parameters for double quantum dot systems.^ 
Chemical potentials of the left and right leads are set to 
be = eV/2 and j-iR = —eV/2, where eV = — j-iR 
denotes the applied bias voltage. 
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FIG. 2: (Color online) Current (a) and differential conduc- 
tance (b) in the parallel (P, solid line) and antiparallel (AP, 
dashed line) magnetic configurations, as well as tunnel mag- 
netoresistance (c), calculated as a function of the bias voltage 
for the parameters: E = 10fc s T, AE = 0, U = 20k B T, 
U = lOOfcsT, p = 0.4, t = 0.25r, T = 5/ieV, a = 0, and 
I = eT/h fa 1.215 nA. 



A. Double dots connected in series, a = 

Let us first consider the situation when a = 0, which 
corresponds to serial geometry of the double quantum 
dot system, see Fig.[TJ In Fig. [2] we show the basic trans- 
port characteristics for the average dot level E = lOfcsT 
and the difference between bare dots' levels AE = 0. In 
the weak coupling regime transport is determined mainly 
by discreteness of the energy dot spectrum and Coulomb 
correlations, which lead to staircase-like current-voltage 
characteristics, see Fig. [21a). For the assumed parame- 
ters, the DQD is empty at low bias and the current is 
blocked below the threshold voltage, irrespective of mag- 
netic configuration of the system. In the blockade regime, 
however, the first-order processes can still contribute to 
the current due to thermal fluctuations. Furthermore, 
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FIG. 3: (Color online) Current (a) and differential conduc- 
tance (b) in the parallel (solid line) and antiparallel (dashed 
line) magnetic configurations, and tunnel magnetoresistance 
(c), calculated as a function of the bias voltage for E — 
— lOfcsT. The other parameters are the same as in Fig. [5] 



in the case of T <C ksT, as considered in this paper, 
the contribution from first-order processes can still be 
larger than that from second-order tunneling (cotunncl- 
ing). Nevertheless, one must bear in mind that in the case 
of deep Coulomb blockade, for instance for E/ksT <C 
and (E + Uo)/ksT 3> 0, the second-order processes be- 
come dominant and must be taken into account to prop- 
erly describe transport properties of the system.— In this 
paper, however, we restrict ourselves to the case when 
transport is mainly governed by sequential tunneling pro- 
cesses. 

When the bias voltage approaches the threshold volt- 
age, the sequential current starts to flow due to onc- 
by-onc electron tunneling through singly-occupied DQD 
states. This leads to the peak in the differential con- 
ductance, see Fig. [Hb). However, when the bias voltage 
increases further, instead of a plateau one observes a drop 
of the current, which leads to negative differential con- 
ductance. This feature appears in both magnetic config- 



urations, see Fig.[2jb). When eV approaches 2E + 2Uq, 
where another electron has possibility to tunnel into the 
DQD system, the current starts increasing further. 

Physical mechanism responsible for the occurrence of 
negative differential conductance follows from the level 
renormalization due to tunneling processes between the 
dots and leads. This renormalization is directly related to 
the real part of the off-diagonal self-energies, see Eq. (|A2|) 
(and also Ref. where the level renormalization in a 
DQD system connected in series and coupled to nonmag- 
netic leads was calculated). Accordingly, the renormal- 
ized level e™ n ' CT of the ith dot for spin a has the following 
form: 



with il^i(x) 



(r^/27r) 
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where 



2nk B T / 

ty(x) is the digamma function. This renormalization lifts 
the initially assumed degeneracy of the two dot's levels. 
The larger separation between these renormalized lev- 
els, the smaller probability of electron tunneling from the 
left to the right dot. We have calculated this separation 
as a function of the bias voltage (not shown) and found 
that the separation increases with increasing bias voltage 
(in the voltage range where negative differential conduc- 
tance appears), and therefore the current decreases with 
increasing bias. After reaching maximum, the level sep- 
aration starts to decrease with a further increase in volt- 
age and negative differential conductance disappears. It 
is also worth noting that, due to the coupling to ferro- 
magnetic leads, the level renormalization becomes spin 
dependent and, consequently, depends on magnetic con- 
figuration of the system, and so does the level spacing. 

The above described level renormalization makes the 
occupation probability of the left dot (QD1) for eV > E 
larger than the occupation probability of the right dot 



(QD2), P\cr.o) > P\o,a) f° r a =t ; 1- Moreover, the rate for 
tunneling between the two dots decreases with increas- 
ing the bias voltage. More specifically, lm(0, cr|/3|cr, 0) de- 
creases, whereas Im(er, 0|/5|0, a) increases, (<r, 0|p]0, cr) = 
((0, cr|/3|(T, 0))*. In the case of tunneling through double 
quantum dots connected in series, the imaginary parts 
of the density matrix, e.g. (0, a\p\a, 0) and (a,0\p\0,a), 
are related to the charge transfer through the system, 
i.e. they are directly related to the current flowJ^ Con- 
sequently, the tunneling of electrons from the left lead to 
the left quantum dot (QD1) is partially suppressed be- 
cause of -P|<r.o> > ^lo.o-) an d due to decreased tunneling 
rates between the two dots with increasing bias voltage 
V in the range \eV\ G (2E, 2E + 2U ). 

The tunnel magnetoresistance as a function of applied 
bias voltage is plotted in Fig.[5{c). Since the conductance 
is larger in the parallel configuration than in the antipar- 
allel one, the corresponding TMR is positive, although 
very small in the voltage range where double occupancy 
is allowed. Moreover, one can note that the sequential 
tunneling TMR is generally smaller than the Jullicre's 
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FIG. 4: (Color online) Renormalization of the dots' levels 
spacing for both spin orientations in the parallel and antipar- 
allel magnetic configurations. The parameters are the same 
as in Fig. [3] 



value of TMR^ TMR Ju11 = 2p 2 /(l - p 2 ) w 0.38 for 
p = 0.4, which is characteristic of single tunnel junc- 
tion or fully coherent transport.— In the case shown in 
Fig.HJc), the TMR reaches local maxima for \eV\ = 2E 
and for \eV\ = 2E + 2Uo, but it is especially enhanced at 
the first Coulomb step, i.e. in the vicinity of \eV\ = 2E. 

In Fig. [3] we present results obtained for DQD con- 
nected in series with average double dot level position 
being negative E = — lOfc^T. The DQD system is then 
singly occupied in equilibrium, and the system is in the 
Coulomb blockade regime as double occupation of DQD 
would cost the correlation energy Uq. Since the levels 
corresponding to E and E + Ug start contributing to cur- 
rent at the same bias (E and E + Uq are symmetric with 
respect to the zero bias voltage), only one step is ob- 
served in the bias dependence of current. Contrary to 
the case presented in Fig. [2j now the current is a mono- 
tonic function of the applied bias voltage, and negative 
differential conductance is not observed, see Fig. [3ja) and 
(b). The difference between the two magnetic configura- 
tions is only slightly visible in the current and differential 
conductance. Accordingly, the TMR above the thresh- 
old voltage is very small, see Fig. [He). At the zero bias, 
however, TMR exhibits a sharp maximum and then drops 
with increasing bias voltage. Furthermore, for voltages 
around the resonance \eV\ ~ 2E a negative TMR is ob- 
served. 

Physical origin of the negative TMR is similar to 
that of negative differential conductance discussed above. 
More specifically, negative TMR follows from the level 
renormalization due to coupling of the dots to external 
leads in the presence of inter-dot Coulomb correlations. 
This level renormalization is spin dependent, so that it 
lifts spin degeneracy of the dot levels. More importantly, 
it is generally different for each dot and therefore modifies 
the renormalized level spacing AE rentCr = e™'' 7 — e^'" ■ 
The renormalization of the level spacing in the system un- 
der consideration is displayed in Fig. [4] For the assumed 
parameters hopping between the dots is like a bottle neck 
for electrons, that controls current flowing through the 
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FIG. 5: (Color online) Tunnel magnetoresistance as a function 
of the bias voltage and the average level position. The other 
parameters as in Fig. [3] 



system. As already stated above, this hopping probabil- 
ity decreases with increasing level separation (for each 
spin orientation). When comparing Fig.[3{c) and Fig.[2J 
one can note that the minimum in (negative) TMR ap- 
pears in the voltage range where the level spacings are 
maximum. Let us consider in more details positive bias, 
eV > 0. In the region where negative TMR is observed, 
the level spacing between left and right dots for the dom- 
inant transport channel (spin-up) in the parallel mag- 
netic configuration is significantly larger than that for 
spin-down electrons in the parallel configuration and also 
significantly larger than the spacing for one of the spin 
channels in the antiparallcl configuration, while it is com- 
parable to the level spacing for second spin channel in 
the antiparallcl configuration. Thus, in the parallel con- 
figuration the spin-down channel, which involves spin- 
minority bands in both leads, takes over control of the 
current, while the dominant spin channel in the antiparal- 
lel configuration involves one spin-majority and one spin- 
minority bands. Consequently, the current is then larger 
in the antiparallcl configuration than in the parallel one, 
which results in negative TMR. 

To support this, we have calculated the relevant occu- 
pation probabilities. We have found that the occupation 
probability of the left dot by a spin-up electron in the par- 
allel configuration is larger than that for the antiparallel 
configuration configuration, P£* > P^-f , whereas oppo- 
site relation holds for spin-down electrons, Pj?, < P£f '■ 
This situation is opposite to that found for the bias re- 
gion where large positive TMR appears. Moreover, the 
occupation probabilities of the right dot by a spin-up or 
spin-down electron in both magnetic configurations are 
small in comparison with those for the left dot. Finally, 
the probability of double occupation of the DQD system 
is rather small, consequently such states are rather irrel- 
evant. 

In Fig. [5] we show the bias voltage and level position 
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FIG. 6: (Color online) The current (a), differential conduc- 
tance (b) in the parallel (solid line) and antiparallel (dashed 
line) configurations, and tunnel magnetoresistance (c) as a 
function of the bias voltage obtained for the parameters: 
a = 1, q = 0.25, E = lOfcgT, while the other parameters 
are as in Fig. [2] 



dependence of the TMR. The position of the average level 
can be changed experimentally by sweeping the gate volt- 
age, so that Fig. [5] effectively shows the bias and gate 
voltage dependence of TMR. One can note that as the 
absolute value of the average level position increases, the 
central maximum in TMR becomes split into two com- 
ponents, while the maximum at zero bias changes into a 
minimum. In the case of negative E, the minimum at 
zero bias occurs roughly when E < — Uq. Furthermore, 
there are also transport regions where TMR changes sign 
and becomes negative. The negative TMR occurs mainly 
for E f=a —Uo/2 and \eV\ « IE. One can also note that 
generally TMR becomes much suppressed for larger bias 
voltages, being close to zero, see also Figs. [He) and^c). 
This implies that the spin polarization of tunneling elec- 
trons is significantly reduced. Such a behavior is opposite 
to that in the case of two strongly coupled dots, where 
TMR in the sequential tunneling regime was found to be 




100 



FIG. 7: (Color online) Current (a) and differential conduc- 
tance (b) in the parallel (solid line) and antiparallel (dashed 
line) magnetic configurations, and tunnel magnetoresistance 
(c), calculated as a function of the bias voltage for a = 0.25. 
The other parameters are the same as in Fig. [6] 



considerable^ 



B. Double dots connected in parallel, a = 1 

In this subsection we consider the case when the dots 
are connected in parallel, a = 1, see Fig. [TJ Virtual tun- 
neling processes between the two quantum dots through 
the leads are now allowed. As mentioned earlier, the 
off-diagonal matrix elements of Tp may be significantly 
reduced due to suppression/cancellation of different con- 
tributions, and hence q < 1 in general. In the following 
we assume q — 0.25. However, we will also analyze how 
the transport characteristics depend on the parameter q. 
For a given q we will examine two cases: symmetric case 
when all dot-lead couplings are the same for given spin, 
i.e., a = 1, and asymmetric case when there is a differ- 
ence in the coupling of a given electrode to the two dots, 
q^0,1. The latter case, referred to as the intermediate 
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geometry, will be analyzed in the next subsection. 

The current, differential conductance and TMR for the 
case when the DQD is empty at equilibrium, E > 0, are 
shown in Fig. [6l The current exhibits a typical staircase- 
like behavior [Fig. [Ufa)], which is also reflected in well- 
defined peaks in the differential conductance located at 
the positions \eV\ = 2E and \eV\ = 2E + 2U , see 
Fig. EJb). The TMR is positive in the whole bias range 
and takes well-defined values corresponding to different 
steps in the current- voltage characteristics, see Fig. EJc). 
In the case of DQD coupled in parallel, the TMR is gener- 
ally larger than in the case of serial connection discussed 
in the previous subsection, especially for large voltages. 
This is due to a different role of interdot hopping in these 
two geometries. As the hopping term is crucial for trans- 
port in serial geometry, it plays a less important role in 
the parallel geometry. Since the hopping parameter is 
independent of spin orientation, it leads to a reduction 
of TMR in the serial geometry in comparison to that in 
the parallel one, especially at large voltages. 



C. Double dots in a general (intermediate) 
geometry 

The situation changes considerably when the dot-lead 
couplings are different, a ^ 1, see Fig. [TJ Figure [7| 
presents the current, differential conductance and TMR 
as a function of bias voltage, calculated for a = 0.25. As 
one can note, the shape of curves describing current and 
differential conductance reveals features obtained above 
for dots connected in series, see Fig. [3J i.e. the nega- 
tive differential conductance and associated non-typical 
Coulomb steps in current. However, the bias depen- 
dence of TMR is now different and is more similar to 
that obtained for parallel geometry, except for voltages 
2E < \eV\ < 2E + 2Uq, i.e. between the resonances, 
where the behavior of TMR is more complex and where 
one finds some oscillations in TMR. Furthermore, for the 
intermediate geometry, the TMR for larger bias voltages 
is quite significant, contrary to the case of serial geom- 
etry, see Figs. [He) and[3tc), and is only a little smaller 
than that in the parallel geometry, see Fig. [5] 

For completeness, in Fig. [5] we show the bias and 
gate voltage dependence of the TMR. One can see that 
TMR displays well defined structures consisting of re- 
gions where it is roughly constant. In these transport 
regions the corresponding current-voltage curves display 
plateaus, whereas at resonances the corresponding TMR 
changes considerably. 

It is interesting to analyze how TMR depends on the 
asymmetry factor a. In Fig. [5] we present the bias de- 
pendence of TMR calculated for different values of a and 
for the case when off-diagonal matrix elements T^ 12 are 
zero (q = 0) and finite (q = 0.25), see Fig. [9ja) and 
(b), respectively. By changing a from to 1, geometry 
of the system continuously changes from serial to paral- 
lel one. Furthermore, with increasing a also the magni- 
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FIG. 8: (Color online) Tunnel magnetoresistance as a function 
of the bias voltage and the average level position, calculated 
for a = 0.25. The other parameters are the same as in Fig. [7] 



tude of the off-diagonal matrix elements is increased, see 
Eqs. © and (J7|). First of all, one can note that TMR 
generally increases with increasing a. This is associated 
with the fact that for a ^ the electrons can tunnel be- 
tween the leads through just a single dot. Accordingly, 
the role of interdot hopping, which is independent of elec- 
tron spin and therefore reduces TMR, is diminished and 
TMR increases. In particular, when crossing over from 
the serial to parallel geometry, the TMR at zero bias in- 
creases, while sharp maxima in TMR at resonance volt- 
ages are transformed into plateaus. This behavior can be 
observed in the case of q = and q = 0.25, see Fig. |5Ja) 
and (b), respectively. There are however some differences 
between these two situations. The virtual tunneling pro- 
cesses between the two dots, described by the nonzero 
parameter q, decrease the TMR at the zero bias and in- 
crease it for large bias voltages, \eV\ > 2E + 2Uo, as 
compared to the case of q = 0. 

Another interesting feature visible in Fig. is that 
TMR for a = 1 at the side plateaus, i.e. for 2E + 2U > 
\eV\ > 2E, has the same magnitude for both q = and 
q = 0.25 cases. This can be clearly shown by deriv- 
ing an approximate analytical formula for TMR in this 
voltage region. At very low temperatures one can ap- 
proximate the Fermi functions by step functions and as- 
sume that the electrons tunnel only from one side to the 
other^ Then, one can show that the TMR for a = 1 
and 2E + 2U a > \eV\ > 2E is given by 



TMR = 



5(1 -P 2 ) 



(12) 



which is exact at zero temperature. From the above for- 
mula follows that TMR in the bias regime under consider- 
ation is independent of the magnitude of the off-diagonal 
matrix elements. This is, however, not true for the cur- 
rent, which for a = 1 is given by 



e r2 
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(13) 




in the parallel and 

/AP = e r (2 + g )(i-p 2 ) 

h 5 + 3p 2 

in the antiparallcl magnetic configurations. 

The interference effects due to virtual first-order 
cesses between the two dots can significantly affect 1 
especially for larger values of q. This is shown in Fii 
which depicts TMR as a function of the bias voltagi 
parameter q, together with various cross-sections. I 
be seen that with increasing the magnitude of th 
diagonal matrix elements, the TMR generally incr 
in the low bias voltage regime, \eV\ < 2E, and for 1 
voltages, \eV\ > 2E + 2U , although the latter di 
dence is not monotonic. However, for bias voltages \ 
transport occurs through charge states with single elec- 
tron on the double dot, 2E + 2U > \eV\ > 2E, TMR 
becomes suppressed with raising q, and for q = 1 we find 
a negative TMR effect. The negative TMR develops ap- 
proximately at the resonance, \eV\ ~ 2E + 2Uq, where 
two-particle states of the double dot system start taking 
part in transport. 

Figure \W\ was calculated for a = 0.25. However, it 
turns our that the spin-dependent transport properties 
may also strongly depend on the geometry of the system. 
This is especially visible in TMR for maximum value of 
the parameter q, q = 1. In Fig. [11] we plot TMR as a 
function of the bias voltage for different values of a and 
for q = 1. As one can note, TMR is rather independent 



FIG. 10: (Color online) (a) Tunnel magnetoresistance as a 
function of the bias voltage and parameter q, calculated for 
for a = 0.25. (b) Cross-sections of TMR for several values 
of q as indicated in the figure. The other parameters are the 
same as in Fig. [6] 



of a at low bias voltages, while for larger bias, \eV\ > 
2E + 2Uq, TMR generally increases with raising a. This 
is however not the case for 2E + 2Uo > \eV\ > 2E, where 
TMR strongly depends on the geometry of the system. 
With crossing over from serial to parallel geometry, the 
magnitude of negative TMR is increased and additional 
maximum develops next to the resonance, \eV\ w 2E + 
2Uq, which transforms into plateau for a close to 1. 

Interestingly, it can be also seen that tunnel magne- 
toresistance is very sensitive to slight changes in the 
asymmetry parameter a, when the latter is close to unity. 
To understand this behavior one needs to realize that the 
bare states of the two quantum dots coupled directly (by 
the hopping term) or indirectly (due to off-diagonal cou- 
pling matrix elements) hybridize in molecular-like states. 
As a result, the bonding and anti-bonding states emerge, 
the widths of which strongly depend on the dot-lead cou- 
pling strengths and geometry of the system. In the case 
considered here, the relative width of the bonding and 
anti-bonding states varies with the parameter a. When 
the difference in the coupling of a given electrode to the 
two dots is reduced, the width of the bonding state in- 
creases whereas that of the anti-bonding state decreases. 
In particular, in the limit of a = 1, the anti-bonding 
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FIG. 11: (Color online) Bias voltage dependence of TMR 
calculated for q = 1 and for different values of the asymmetry 
parameter a, as indicated in the figure. The other parameters 
are the same as in Fig. [2] 
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FIG. 12: Graphical representation of the self-energy SL' { L' , 
(a). The summations are over lead and spin degrees of free- 
dom a = L,R and a' =Ti4- Each tunneling line carries lead 
index a, spin a and frequency u>. 



negative differential conductance vanish. The above ef- 
fects may be restored in an intermediate geometry and 
strongly depend on the magnitude of the virtual pro- 
cesses between the two dots. Furthermore, in the case 
when virtual processes are maximal, we have found a 
strong dependence of the TMR on the geometry of the 
system, especially for geometries very close to the parallel 
one. 



state becomes totally decoupled from the leads, while 
the bonding state acquires width of the order of 2T. 
In this limit the above mentioned features of the TMR 
disappear and tunnel magnetoresistance is constant for 
2E + 2U > \eV\ > 2E. In other words, the high sensitiv- 
ity of the TMR with respect to the system's geometry is 
associated with the fact that the interference conditions 
for electron waves transmitted through the two dots be- 
come modified when changing a. 



IV. SUMMARY AND CONCLUSIONS 

We have analyzed the spin-polarized transport proper- 
ties of double quantum dots weakly coupled to each other 
and to external leads. Using the real-time diagrammatic 
technique we have calculated the conductance and tun- 
nel magnetoresistance in the parallel, serial and interme- 
diate geometries of double quantum dots. Moreover, we 
have taken into account the effects of virtual tunneling 
processes between the two dots taking place through the 
states in the leads. Such processes are absent in serial 
geometry and become maximum for parallel geometry. 

In the case of double quantum dots coupled in scries 
we have found a negative tunnel magnetoresistance at 
the resonance and negative differential conductance for 
transport voltages where single-particle double-dot states 
take part in transport. On the other hand, for paral- 
lel geometry of the system, both the negative TMR and 
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APPENDIX A: EXAMPLES OF THE 
FIRST-ORDER SELF ENERGIES 

In order to analyze the transport properties one needs 
to calculate the respective self-energies using correspond- 
ing diagrammatic rulesj 23 i 34 In the sequential tunneling 
regime only the first-order self-energies determine the 
transport characteristics. Here, we present explicitly two 
examples of first-order self-energies. Generally, the self- 
energies are complex - their imaginary part may be re- 
lated to transition rates, whereas the real part may be 
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associated with with various renormalization effects. The displayed in Fig. 1121 while analytically it is given by 
graphical representation of the self-energy e!^' !!^' ! is 
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FIG. 13: Graphical equation for the self-energy Zgy 



with 



P. 



±,<7, 



x) = ±- 



2tt 



Re* 



l 2nk B T 



- log 



2nk B T 
(A3) 



In the above expression the states |2, 0), |0, 2) denote dou- 
bly occupied first and second dot, respectively, and, be- 
cause of large intradot Coulomb repulsion, are only con- 
sidered as virtual ones (intermediate states). In Fig. [TBI 
we also show an example of self-energy including off- 
diagonal matrix elements, Fg, YuZ t'oj, it is given by 

E |K!fe?> = E KT(^) - 

+ Pif (e 2ff + C/o) + nr 7 +' 2 CT (e 2 a + C/o) } , (A2) 



7^(*) = ^/ ± (*-A*a) 1 (A4) 



where f + (x) stands for Fermi-Dirac distribution, 
f~(x) = 1 — f + (x), and *(x) denotes the digamma func- 
tion. Here, D is the cut-off parameter which can be iden- 
tified with on-level Coulomb repulsion U. It is worth to 
note that the self-energies do not depend on the cut-off 
parameter because the terms depending on D cancel in 
pairs, see Eqs. ([AT]) and (|A"2|) . 
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